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ABSTRACT 

We address the origin of the 'downsizing' of elliptical galaxies, according to which 
the stars in more massive galaxies formed earlier and over a shorter period than those 
in less massive galaxies. We show that this could be the natural result of a shutdown 
of star formation in dark matter haloes above a critical mass of ~ 10 12 M Q . This is 
demonstrated using a semi-analytic simulation of galaxy formation within the standard 
hierarchical scenario of structure formation. The assumed threshold mass is motivated 
by the prediction of stable shock heating above this mass and the finding that such a 
shutdown reproduces the observed distribution of galaxies in luminosity and colour. 
The shutdown at a critical halo mass introduces a characteristic stellar mass for the 
transition of galaxies into the 'red sequence' of the galaxy colour-magnitude diagram. 
Central galaxies of haloes that are more massive today have reached this mass earlier 
and can therefore grow further along the red sequence by dry mergers, ending up more 
massive and containing older stars. Small galaxies formed in haloes below the critical 
mass can shutdown late, when they fall into haloes above the critical mass and become 
satellites. While our semi-analytic simulation that incorporates an explicit shutdown 
reproduces downsizing as inferred from the stellar ages of ellipticals, we explain why 
it is much harder to detect downsizing using the mass functions of different galaxy 
types. 
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1 INTRODUCTION 



2005] : llleda et all 120031; iHeckman et alj 12004 iBarger et ail 
2005t iHasinger et alj|2005l ) . 



The term 'downsizing' was coined by ICowie et alj (|l996l ) 
to describe the decline with time of the Tf-band rest- 
frame luminosity of the galaxies with the highest spe- 
cific star formation rate M s tar/M s tar (SFR) as ob- 
served in the redshift interval 0.2 < z < 1.7. 
This 'do wnsizing in time' ha s been confirmed by later 



studies dGuzman et al 
Kodama et all 12004; 



1997 



Noeske et all 



20071) . 



Brinchma nn fc Ellis 
Juneau et alj 120051 : | Bell et alj 



2000 



2005 



It can be viewed as part of a more 
general phenomenon valid robustly across the Hubble se- 
quence of galaxy types, in the sense that galaxies of later 
type, which are typically less massive, are known to have 
formed their st ars most efficiently at later times and over 
longe r periods (|Searle et alj 1 19731 : iTinslevI 1 19731 ; ISandagej 
1 19861 ) . A possibly related downsizing is evide nt in the ac- 
cretion histories of supermassive black holes (|Steffen et al.1 



Another form of downsizing , termed 'archaeological 
downsizing' (I Thomas et alj [2005), is inferred from the 
stellar populations of today's galaxies. Star formation 
histories derived from observed line indices and abun- 
dance ratios using stellar evolution models reveal a strong 
correlation between mean stellar age and galacti c stel- 
lar mass both in elliptical galaxies dNelan et al.1 120051 ; 
iThomas et~aT1 120051 ; iGraves et ail 120071 ) and in the gen- 
eral population of galaxies from the large Sloan Digital 
Sky Survey (SPSS ; iHeavens et al]|2004 I Jimenez et alll2005l ; 
iPanter et alj|2006l ) . In the current study we focus on archae- 
ological downsizing, as we address red galaxies that show 
only little active star formation at present. 

At a first glance, the observed downsizing seems to be in 
conflict with the standard wisdom of hierarchical structure 
formation, according to which smaller objects collapse ear- 
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lier and gradually assemble into more massive objects. This 
apparent conflict is relaxed when realizing that star forma- 
tion and gravitational assembly are two distinct processes. 
In particular, stars can form first in the small building blocks 
of today's massive galaxies without violating the assembly 
hierarchy of each galaxy. In fact, if gas processes limit galaxy 
formation to dark matter haloes above a minimum mass, a 
certain downsizing trend arises naturally fr om the hierar- 
chica l dark matter assembly process itself (|Neistein et al.l 
l2006h . Still, the origin of downsizing as observed requires a 
more quantitative theoretical understanding. 

Recent dev elopments in t he modelling of galax> 
formation (e .j 



slopments m t ne m odelling ol galaxy 
Cattaneo et~ai] 120061 ; iBower et all 120061 : 



ICroton etT al, 2006) have been driven by the realization that 
galaxies are divided into two major distinct types: blue 
star-forming late-type galaxies in low-density environments, 



dLin et al. 19971. 


19991; 


Im et al.l 2001 ; IStrateva et al. 200 ll 


Kauffmann et al. 


2003 


; Baldrv et al. 2004; Balogh et al.l 


2004 Bell et al. 20041: Hoe 


z et al.ll2004l; IWeiner et alj|2005l; 


see Dekel & Birnboimll2006 


for a summary of the bimodality 



phenomenon) . Combined with the fact that blue galaxies are 



confined below a critical stellar mass of MS 



3 x 1O 1O M 

while the population above this mass is dominated by red 
galaxies, the colour bimodality itself is evidence for down- 
sizing. It indicates that the most massive galaxies, the red 
ones, have converted their gas into stars several billion years 
ago, while less massive galaxies, which are mostly blue, are 
still making stars. 

Downsizing is not only a reflection of the galaxy bi- 
modality. It is also valid s eparately within each of the two 
popu l ations, the red one (|Nelan et al.ll2005l ; iThomas et al 



2005 



2006 



Graves et all 120071) and the blue one (|Drorv et al 



Noeske et all 120071 ). Less massive galaxies contain 



younger stars even within the same spectral type. Neverthe- 
less, we are tempted to consider the possibility that down- 
sizing and bimodality originate from the same underlying 
physics. 

The existence of a characteristic mass for mas- 
sive galaxies emerges from the the competition be- 
tween the gravitati onal dynamical time a nd the radia- 
tive cooli n g time ([flees fc Ostrikerl 1 19771; iBinne yl 1 19771 ; 
Silkl Il977l; IWhite fc Reesl Il978l ; iBlumenthal et all Il984h . 
Birnboim fc Dekell (|2003h and iDekel fc Birnboiml (|2006l ) 



showed that a stable shock can expand to the virial ra- 
dius and be supported against gravity only in haloes above 
a critical mass of M s h oc k ~ 10 12 Mq, producing a hot 
medium at the virial temperature of T ^ 10 6 K. In smaller 
haloes, efficient cooling does not allow an extended sta- 
ble shock, so the accreting gas flows cold into the cen- 
tre, leading to the buildup of discs and efficient star for- 
mation. This phenomen on was detected i n parallel in cos- 
molo g ical simulations dKeres et al.l 1200 5l ; iBirnboim et al.l 
120071 ; ICattaneo etTafl 120071 ). confirming the existence of 
a threshold mass that is quite independent of redshift. 
IDekel fc Birnboiml (|2006l ) argued that the transition from 
cold flows to the shutdown of gas supply and star forma- 
tion at JV/shock, which is predicted to be more pronounced 
after z ~ 2 — 3, naturally leads to many of the observed 
features associated with the galaxy bimodality. In partic- 
ular, the bimodality scale M££a is the typical stellar mass 
of the central galaxy in a halo with mass M s hock- The im- 



plications of such a sc enario have been studi e d using semi- 
analyt ic si mulations bvlCattaneo et al.l (|2006l ). ICroton et al.l 
(|2006l ) and lBower et all (|2006l ). each applying a slightly dif- 
ferent shutdown procedure driven by halo mass (plus, some- 
times, bulge mass or central black hole mass) . The formation 
of the red sequence through the quenchi ng and reddening 
of b lue galaxies has als o been studied bv lBell et al.l (|2004 ) 
and iFaber et al.l (|2007l ). Here we highlight a particular as- 
pect of the quenching mechanism, associated with the fact 
that the critical mass for quenching is roughly constant after 
z ~ 2 - 3. 

In this article we demonstrate that the shutdown of star 
formation above a critical halo mass M cr i t ~ M s hock at z < 3 
accounts for downsizing in the red galaxy population. The 
key idea is that central galaxies stop making stars and enter 
the red sequence at a critical stellar mass M s ™5 of the or- 
der of the baryonic mass in a halo of mass M s hock- Galaxies 
with higher final stellar mass reach M s c t "' earlier and there- 
fore have more time to grow through dissipationless ('dry') 
mergers along the red sequence. 

There are three complications to this simple picture, 
which we discuss in this article. First, satellite galaxies shut 
down earlier than central galaxies with the same mass, and 
they often do it before they have reached Mg t "* (e.g. because 
their halo has merged into another one with mass > -Mshock). 
Second, the expected penetration of narrow cold streams 
through the hot media of massive high-redshift haloes ele- 
vates the effective M cr it to values above M s h o c k at z ^ 3 
(|Dekel fc Birnboiml |2006| ; ICattaneo et al.ll2006l ). Third, ad- 
ditional shutdown mechanisms may be at work, e.g. the ex- 
haustion of gas after major mergers and/or feedback from 
the central black hole. We have attempted to mimick these 
effects by intoducing a further shutdown criterion related 
to the bulge-to-disc ratio. The inclusion of a bulge-to-disc 
ratio criterion as a sufficient condition for quenching in ad- 
dition to the halo- mass criterion improves the fit to the joint 
colour-magnitude distribution of galaxies by accounting for 
the presence of red galaxies in haloes below M cr j t , but its 
overall effect is fairly minor (|Cattaneo et alj|2006l . Fig. 9). 

Downsizing has also been involved to explain the way 
the stellar-mass functions of galaxies of different colors or 
morphological types evolve with time. In particular, while 
the comoving number density of blue galaxies has not 
changed much since z ~ 1, the number of red galaxies 



than at the 


bright end dRudnick et al.ll2003l; Bell et alJ 


2004 


Drorv et al. 


2004 Bundv et al. 20051; iDrorv et al.l 


2005 


Borch et al. 


120061; iPannella et all 120061; IFaber et al.l 120071) 



We argue below that the interpretation of this behavior as 
downsizing is hampered by severe measurement uncertain- 
ties (see the Appendix) and conceptual misunderstandings. 

The paper outline is as follows. In § 2, we describe our 
model of galaxy formation explored by means of a semi- 
analytic modelling technique applied to m erger trees from 
a cosm ological N-body simulation, following ICattaneo et al.l 
(|2006l ). In § 3, we elaborate on how the shutdown of star 
formation leads to downsizing. In § 4, we show how the 
shutdown produces downsizing in the buildup of the red se- 
quence, in the sense that the high-mass end is populated 
at earlier times. In § 5, we perform a comparison with the 
archeological do wnsizing inferred from the stellar ages in el- 
liptical galaxies (|Thomas et al.l [20051 ) . In § 6, we consider 
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the characteristic times of galaxy formation in terms of star 
formation history and mass assembly history, and use this 
analysis to address the apparent conflict between downsizing 
and the hierarchical model. In § 7, we explain how downsiz- 
ing in star formation and upsizing in mass assembly combine 
to produce the evolution of the mass functions of red and 
blue galaxies, and why it is hard to detect downsizing us- 
ing the mass functions of different galaxy types. In § 8, we 
summarize the conclusions of the article. The downsizing in 
blue galaxies will be addressed in a separate publication. 



2 THE GALAXY FORMATION MODEL 

GallCS (Galaxi es In Cosmological Simulations; 
lHatton et afll2003t ) is a method to simulate the formation of 
galaxies in a ACDM Universe. It combines high-resolution 
cosmological N-body simulations of the gravitational 
clustering of the dark matter with a semi-analytic (SAM) 
approach to the physics of the baryons (gas accretion, 
galaxy mergers, star formation and feedback). The version 
of GallCS used here, as well as the adopte d values of the 
param eters, are the same as described in ICattaneo et alj 
l)2006l . and references therein). We summarize here the 
relevant issues. 



2.1 Dark-matter simulation 

The cosmological N-body simulation that follows the hierar- 
chical clustering of the dark-matter component has been car- 
ried out with the parallel tree code. The assumed cosmologi- 
cal model is a flat ACDM Universe with a cosmological con- 
stant of f^A = 0.667, a Hubble constant of Ho — 66.7 km s , 
and a ACDM power spectrum of initial fluctuations normal- 
ized to erg = 0.88. The computational volume is a cube of 
side (150 Mpc) 3 with 256 3 particles of 8.3 x 1O 9 M each 
and a smoothing length of 29.3 kpc. The simulation pro- 
duced 100 output snapshots spaced logarithmically in the 
expansion factor (1 + z) _1 from z = 35.59 to z = 0. 

We ha ve analysed each sn apshot with a friend-of-friend 
algorithm (|Davis et al.l [1983) to identify virialized haloes 
containing more than 20 particles. The minimum halo mass 
is thus 1.65 x lO n M0. The three global properties charac- 
terizing each halo to be used in the SAM are the virial mass, 
the virial density, and the spin parameter. Merger trees are 
constructed by linking the haloes identified in each snapshot 
with their progenitors in the previous snapshot, that is, all 
predecessors from which the halo has inherited one or more 
particles. We do not use substructure information from the 
N-body simulation. Instead, once a halo becomes a subhalo 
of another halo, we switch from following its evolution with 
the N-body integrator to an approximate treatment based 
on semi-analytic prescriptions. 



2.2 Semi-analytic modeling of gas processes 

A newly identified halo is assigned a gas mass based on a 
universal baryon fraction f2b/£2o = 0.135. The conditions 
for the formation of a galaxy at its centre are that the halo 
is gravitationally bound and that its angular momentum 
parameter is A < 0.5. We have simulated a 'standard' and 



a 'new' model, which differ in the way gas from the halo 
accretes onto the central galaxy. 

In the 'standard' model, the initial gas distribution 
is a singular isothermal sphere, at the virial temperature, 
truncated at the virial radius. The cooling time of the hot 
gas is computed usi n g the radiative cooling function of 
[Sutherland fc Dopital (|1993h . The gas for which both the 
cooling time and the free fall time are shorter than the 
timestep At between two N-body snapshots is accreted onto 
the central galaxy during that time interval. The transfer of 
baryons from the halo to the galaxy is accompanied by an 
inflow of hot gas to keep the hot gas distribution a truncated 
singular isothermal sphere. The numerical accuracy and the 
self-consistency of this cooling algorithm have been verified 
by comparison with a cosmological hydrodynamic simula- 
tion based on the same physics and the same dark-matter 
realization (|Cattaneo et alj|2007f ). 

In the 'new' model, we introduce the proposed 
quenching above a th reshold halo mass, following 
iDekel fc Birnboiml <|2006h . Once a halo grows above 
Afshock ~ 10 12 Mq, the halo gas is assumed to be shock- 
heated to the virial temperature. Unlike previous models, 
once heated, the hot gas is kept hot forever. This can be due 
to several different long-term quenching mec hanisms, such 
as se l f-regulated radio-AGN fee dback (e.g. ICroton et al.l 
120061 ; ICattaneo fc T cvssicr 2007, a nd re ferences therein), 



shocked accretion ( Birnboim et al. | 120071). or gravitational 
quenching by clumpy acc retion ( Dekel fe Birnboiml 120071 ; 
iKhochfar fc Ostrikerll2007l ) . The critical mass for quenching 
is assumed to be 



Merit = M shock x max{l, 10 



1.3(z-z c ) 



} • 



(1) 



where M s hock ~ 1O 12 M0 is the shock-heating scale and the 
term oc 10 1,3z accounts for the pe netration of cold streams 
in massive haloes at high redshif t (|Dekel fc Birnboiml 120061 . 
Fig. 7). In lCattaneo et all (|2006l ). we set the parameters to 
-M s hock = 2x 10 12 Mq and z c = 3.2 by optimizing the way the 
new model fits the colour-magnitude distribution at z ~ 
and the Lyman-break galaxy luminosity function at z ~ 3. 

In the standard model, there is no explicit shutdown of 
cooling, and the mass and luminosity of the central galaxy 
grow roughly linearly with halo mass over the entire range 
resolved by the N-body simulation. Shutting down the hot 
accretion mode when Mhaio > M cr i t introduces a character- 
istic mass (Fig.[TJ and luminosity, which mark the separation 
between blue and red galaxies. 

Initially, all galaxies are assumed to form at a distance 
n — from the centre of their halo. When two or more 
haloes merge, their galaxies are repositioned at a distance 
77 from the centre of the new halo, determined such that 
r f — > n for Mf/Mi — » 1 and 77 — » r v i r for Mf/Mi — » 00 
(where Mi and Mf are the halo mass before and after the 
merger and r v i r is the virial radius of the merged halo. Halo 
mergers rapidly create haloes with more than one galaxy, 
which we identify as galaxy groups and clusters. However, 
only the central galaxy is allowed to accrete gas from the 
halo (prior to quenching). 



2.3 Morphologies 

GallCS models a galaxy with three components: a disc, a 
bulge and a transitional 'starburst' component. Each com- 
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New Model 
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Figure 1. Galaxy stellar mass versus halo mass in the standard model (left) and the new model with shutdown (right). In the new 
model, the general growth of M s t ar with M^aio breaks off at M^aio ~ A^ock ~2x 10 12 Mq. Differences at Afjjalo < Af s hock are due to 
the additional shutdown of accretion onto the galaxies that have become bulge-dominated. Objects with high iWhai but low M gal are 
mostly satellites that have fallen into larger haloes. 



ponent may contain stars and cold gas, while the hot gas 
is treated as a component of the halo. Only the disc of the 
central galaxy accretes gas from the halo. The bulge grows 
by mergers and by disc instabilities. While stars are tras- 
ferred directly from the disc to the bulge, the gas in transit 
from the disc to the bulge is assumed to pass through an 
intermediate starburst component, where its star formation 
timescale decreases by a factor of ten with respect to that 
of the bulge (Section 2.4). Stars formed in the starburst 
are moved to the bulge after they have reached an age of 
100 Myr. The only gas in the bulge is that from stellar mass 
loss (§ 2.4). 

The disc profile is assumed to be exponential, with its 
radius determined by conserv ation of angular momentum. 
The bulge is assumed to have a iHernauistl (|l99Ch profile and 
its radius is determined based on an energy conservation 
argument. 

We now elaborate on the two mechanisms for trans- 
ferring baryons from the disc to the spheroid. First, disc 
instabilities transfer matter from the disc to the spheroid 
(stars to the bulge, gas to the starburst) until the bulge 
is massive enough to stabilize the disc. The stability cri- 
terion is w ro t < 0.7u c , where D rot and w c are the disc's ro- 
tation and circular v elocity at the disc's half mass radius 
(Ivan den Boschlll998l). Sec ond, dynamical friction, modelled 



M star = — — (1 + z) 

P* £(lyn 



(2) 



Hattonetain fa. drives galaxies to the centre of 



their dark matter halo and is the dominant cause of galaxy 
merging, although we also include satellite-satellite encoun- 
ters. The fraction of the disc that is transferred to the bulge 
in a merger grows with the mass ratio of the merging galax- 
ies. It ranges from zero for a very minor merger to unity 
for an equal mass merger. This simplified picture of the 
dynamics of morphological transformations provides results 
consistent with the key observational constraints such as 
the Faber- Jackson relation an d the Fundamental Plane of 
spheroids ijHatton et al.ll2003T ). 



2.4 Star formation and feedback 

The star formation law is the same for all components: 



The mass of cold gas, M 00 ld, refers to the component in ques- 
tion, and tdyn is the dynamical time (corresponding to half a 
rotation for discs and half a crossing time for bulges) . In the 
standard GallCS model, star formation is activated when 
the gas surface density > 20 m p cm (m p is the pro- 

ton mass) . Furthermore, by imposing a minimum halo mass 
of Mhaio = 1-65 x 1O 11 M0 due to the N-body resolution, we 
effectively assume that there is no star formation in haloes 
below this mass, potentially mimicking the effects of more 
aggressive supernova feedback in small haloes (see below). 
Moreover, in the new model we assume a complete shutdown 
of star formation in haloes with mass > M CT n by removing 
all the cold present in any component of these galaxies and 
by moving this gas to the halo's hot component. 

The sta r formation efficiency p arameter is assumed to 
be /3» = 50 (jGuiderdoni et al"1ll99ct ). the same for all com- 
ponents. We have a* = in the standard GallCS model but 
allow an enhanced star-formation rate at high redshift by 
adopting a* — 0.6 in the new model, to mimic the rapid ac- 
cretion by cold flows. This improve s the fit to the lumino sity 
function of Lyman-break galaxies (jCattaneo et al.ll2006l) . 

Gas in transit from the disc to the bulge, whether 
by mergers or disc instabilities, passes through a starburst 
phase in which the SFR grows by a factor of 10 based on the 
dynamical time of the bulge. This high SFR is obtained by 
assuming that the starburst radius is ten times smaller than 
the final bulge radiu s (see the hydrodynamic simulation in 
ICattaneo et al.|[2005l ). As we anticipated when we described 
the different components of a galaxy, the stars formed in 
the starburst component stay in the starburst component 
for 100 Myr and then are moved to the bulge. The proper- 
ties of low redshift galaxies are insensitive to the starburst 
star formation timescale as long as it is much shorter than 
the Hubble time. 

We assume a iKennicutti (|l983T l initial mass function. 
Stars are evolved between snapshots using substeps of at 
most 1 Myr. During each sub-step, stars release mass and 
energy into the interstellar medium. Most of the mass comes 
from the red giant and the asymptotic giant branches of 
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Standard Model No Additional Bulge Mass Criterion New Model 
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Figure 2. The the u — r colour distribution in r-band magnitude bins for the SDSS data jBaldry et al I I2004 black smoothed histogram) 
and for three different versions of our SAM (red curves). Left: the standard model without an explicit shutdown of central galaxies but 
with a standard shutdown of gas accretion onto satellites. Middle: A model with shutdown of st ar formation once M^a io > Merit' Right: 
Same as in the middle, with the additional shutdown when Mb u i ge > M s tar/2 (reproduced from lCattaneo et al,|[2006i) . 



stellar evolution, while most of the energy comes from shocks 
due to supernova explosions. The enriched material released 
in the late stages of stellar evolution is mixed with the cold 
phase, while the energy released from supernovae is used to 
reheat the cold gas and to return it to the hot phase in the 
halo. Reheated gas is ejected from the halo if the potential 
is shallow enough. The rate of mass loss through supernova- 
driven winds M w is determined by the equation 

1-2 

-MwlW — eSN^SN-EsN-Mstar, (3) 

where _Esn = 10 erg is the energy of a supernova, t]sn = 
0.0093 is the number of superno vae for 1 Mq of sta rs formed 
and Vesc is the escape velocity l|Dekel fe Silk|[l98r3 ). In Gal- 
ICS, we use iw — 1.84u c for discs and v csc — 2a for bulges 
and starbursts. The supernova efficie ncy 6sn — 0-2 is similar 
to those commonly ado pted in SAMs (|Somerville fc Primackl 
ll999l ; ICole et alj|2000h . 



2.5 Stardust 

The STARDUST model (|Devriendt et al.lll999i ) provides the 
spectrum of a stellar population as a function of age and 
metallicity and includes a phenomenological treatment of 
the reprocessing of stellar light by dust. GallCS uses STAR- 
DUST to compute the spectrum of each component and to 
output galaxy magnitudes. 

Dust absorption is computed with a phenomenological 
extinction law that depends on the column density of neu- 
tral hydrogen, the metallicity of the obscuring material and 
the randomly selected inclination angle (for spirals) . The re- 
emitted spectrum is the sum of four templates (big and small 
carbon grains, silicates, and polycyclic aromatic hydrocar- 
bons). Their relative weights are chosen to reproduce the 



relation between bolometric luminosity and infrared colours 
observed locally in IRAS galaxies. In GallCS we only con- 
sider each component's self-absorption, e.g. the dust in a 
starburst only obscures the starburst's light, not the old 
bulge stellar population. 



3 SHUTDOWN AND DOWNSIZING 

In the new model, the accretion of gas onto the central 
galaxy is shut down once Mhaio > Af cr it- Both in the stan- 
dard and the new model, as common in other SAMs, the 
accretion of gas onto satellite galaxies is not allowed^ Fig- 
ure [1] shows that the stellar mass keeps pace with the halo 
mass up to M cr i t and then breaks off in the new model. It 
also shows the effect of preventing the accretion of gas onto 
satellite galaxies, which are shown by the point clouds in the 
lower right part of the diagrams. 

In order to maximize the quenching effect, whenever 
the shutdown criteria are satisfied we actually shut down 
star formation altogether, that is, in addition to having the 
halo gas shock heated and to keeping it hot, we assume that 
the cold gas is removed from the galaxies and that it is added 
to the hot halo. This can be achieved in cent ral galaxies by 
thermal evaporation (|Nipoti fc Bin ncv 2007j) , and in satel- 
lite galaxies by ram-pressure str i pping (e.g. iBahcalll Il977l ; 
lGallagherlll978h . ICattaneo et al.1 (|2006r ) demonstrated that 
this robust shutdown leads to a good fit to the observed 

1 Admittedly, this procedure may be too strict and should be 
revised, especially in small haloes with no hot medium, where 
anyway it may be hard to distin guish between central galaxies 
and s atellites (e.g. ICattaneo et al,ll2007l ; also sec McCarth y et al] 
l200Tt) . 
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Figure 3. Colour-mass diagram for the new model with shutdown. The 'green valley' separating the red sequence from the blue cloud 
is marked by the black line. Blue symbols mark galaxies in which the SFR has not been explicitly shut down. Red symbols denote 
central galaxies in which the SFR has been shut down because their halo has grown above the critical mass. Orange symbols show the 
population of satellite galaxies in haloes above the critical mass. Green symbols are central galaxies of haloes with M^aio < A^criti where 
the formation of stars has been shut down because the bulge has become the dominant component. Black symbols refer to galaxies that 
do not accrete gas simply because they are satellites. In this diagram galaxies have been coloured according to their status at z = 0, and 
not according to the first mechanism that has intervened and caused the shutdown of star formation in the history of a given galaxy. 
The red parallelograms show the four bins along the red sequence from which we have selected the galaxies for Figs. [4] and [5] 



colour- magnitude distribution (Fig. [2j centre). In the ab- 
sence of such an explicit shutdown the standard model pre- 
dicts too many luminous and blue galaxies (Fig. left). 

The only slight disagreement between the predictions 
of the simple model of shutdown by halo mass (red curves 
in the panel marked 'no additional bulge mass crite- 
rion') and the SDSS data (black smoothed histograms) 
is a small excess of blue galaxies at —22.25 ^ M r , which 
arises from the absence of quenched central galaxies with 
M s t a r ~ 1 — 2 x 10 11 , since their haloes are below the 
shock-heating scale. This is the range where the ellip- 
tical population is dominated by rather discy configura- 
tions, which do not show X-ray e mission in excess o f 
the contribution of discrete sources (|Bender et al.l [T989T ) . 
These galaxies could have been possibly quenched by 
another mech anism, such a s quasar feedback after gas - 
rich mergers jSpringel et al. | l2005al ; IHopkins et ail l2007h . 
following the notio n that di s cy el lipticals emerge from 
wet mergers (e. g. ICox et al.l 12006 ) , wh i ch also trigger 
quasar activity (IToomre fc Toomre 19721; Cattaneo et al.l 



Il999l ; iKauffmann fc Haehneltll2000l ; ICattaneo et al J l2005l a.b 



ISpringeT et al.l l2005bl ; IHopkins et all 120051 ). We do not at- 
tempt here to model quasar feedback in detail. We simply 
mimic its effect by stopping the accretion onto a galaxy 
whenever it is dominated by its bulge component, even if 
it is still in a halo with Mhaio < Merit- This can occur only 
after a major merger or a sequence of minor mergers since 
the disc instability criterion does not permit the growth of 
bulges with Aibuigc > M s tar/2 (see Section 2.3). This ad- 
ditional quenching criterion slightly improves the fit to the 
SDSS colour-magnitude distribution near M r ~ —21.5 (the 
'new model' in Fig[2J ■ Note that the cold gas is not removed 
when a galaxy is quenched by this criterion alone. While 
this is clearly a simplified ad hoc implementation of quasar 
feedback, it has negligible effects at the bright end and at 
the faint end, and it hardly affects the main results of our 
current analysis. 

One could conceive an alternative scenario where quasar 
feedback af ter galaxy mergers i s the main shutdown mech- 
anism (e.g. IHopkins et a.1.1 120071 ') . although it would have a 
hard time explaining the long-term maintenance of quench- 
ing. There is some observational evidence for such a mecha- 
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Figure 4. Growth of halo mass (top) and stellar mass (bottom) for 24 galaxies randomly sampled from the four mass bins indicated 
(see Fig. [3}. The growth curves refer to the most massive progenitors. The bold black line in the upper panels is M cr ;t(z) above which 
cooling and star formation are shut down (Eq. 1). After z ~ 3, M cr ; t ~ M Bllock ~ 2 X 10 12 Mq. After M ha i D has crossed this threshold, 
the central galaxy ceases to form stars and passively turns red and dead. Growth of the stellar mass after this point occurs only via 
mergers. The sudden jumps in Afhalo occur when a galaxy becomes a satellite in a larger halo. The passive fading of giant ellipticals 
after z ~ 3 — 4 is well described by the simple relation Mb ~ — 21.5 — z ±0.5 (solid line between the two dotted lines corresponding to 
±0.5 mag) 
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nism in action (jSchawinski et afll2007h . Modelling this spe- 
cific scenario is beyond the scope of the present paper, where 
we focus on the effect of halo quenching on downsizing. 

In summary, the new model stops the accretion of gas 
onto a galaxy when at least one of the following three con- 
ditions is satisfied: (a) Mhaio > M cr it, given by Eq. 1, or 
(b) Mbuige > M s tar/2, (c) the galaxy has become a satel- 
lite. When the quenching is by condition (a), the cold gas is 
explicitly removed, but even in the other cases the starved 
galaxies rapidly exhaust their gas and stop making stars. 
Fig. [3] illustrates via different colours the role of these three 
different mechanisms in producing the galaxy bimodality. 
We see that halo-mass quenching (a) is the main mechanism 
for M s tar ~2x 10 11 M o , and that the bulge quenching crite- 
rion (b) plays a role near M ata r ~ 10 n Mo, while most low 
mass galaxies are both quenched by halo mass (a) and satel- 
lites (c). Only 3% of the red sequence at Af Bta r > 10 10 M Q 
is made of galaxies that are quenched only because they are 
satellites (c) . Although 57% of the red sequence is composed 
of galaxies with Mbulge > M star /2, only 24% of red galaxies 
are quenched by criterion (b) only; ~ 10% of the galax- 
ies quenched by criterion (b) only are red because they are 
dusty. These galaxies are remnants of recent gas-rich merg- 
ers. They have ceased to accrete gas, but they have not yet 
ceased to make stars. They account for ~ 2% of the red se- 
quence at M star > 1O 1O M0. About 8% of the red sequence 
at Mstar > 10 10 Mo is made of galaxies that have not been 
explicitly quenched by any of the three mechanisms above 
(blue points in Fig. [3J . The breakdown of this 8% is as fol- 
lows. About 7% are spiral galaxies that have naturally run 
out of gas without an obvious external trigger, and the other 
1% are red because they are dusty. In conclusion, 3% of the 
red sequence, a rather low fraction, is red because of dust. 

The shutdown of star formation when the halo mass 
grows above M cr it introduces the notion of an 'entry mass' 
to the red sequence, M^', that is t he typical stellar m ass of 
the central galaxy of an M cr it halo (|Faber et al.ll2007l l The 
entry mass is predict ed to be constant after z ~ 2 and higher 
at higher redshifts (|Dekel fc Birnboimll2006t) . Its value at 
2 2 is under 10% of Mshock ~ 10 12 Mq, given the univer- 
sal baryonic fraction of ~ 0.15 and assuming that at least 
half of the gas is prevented from making stars due to stellar 
feedback. 

We illustrate the connection between shutdown, entry 
mass and downsizing by comparing the halo and stellar mass 
growth histories of galaxies with different final stellar masses 
(Fig. [4| . We consider four mass bins along the red sequence 
in the colour-mass diagram (corresponding to the four red 
parallelograms in Fig. [3} and select six galaxies at random 
from each bin. The mass intervals corresponding to the four 
bins cover a factor of forty in stellar mass, from ~ 0.2M s c t r a ' 
to ~ 8M™r- 

All galaxies display an initial phase, in which the stel- 
lar mass grows rapidly, followed by shutdown and a passive 
phase, in which the stellar mass remains essentially constant 
(Fig. [4| . The transition between the two phases is accompa- 
nied by a very fast reddening of the U — B galaxy colour. This 
reddening is recognizable as a vertical jump upward in a di- 
agram of colour vs. stellar mass (Fig. [5j left column). It can 
also be seen in the galaxies' evolutionary colour-magnitude 
tracks (Fig. [5] right column), as it corresponds to the point 
where the galaxies start fading and reddening. In 23 out of 



24 galaxies, the transition from blue to red coincides with 
the time when the halo mass crosses M cr i t . A slight decrease 
of the stellar mass in the passive phase is due to the death 
of old stars not being replaced by new ones. Growth in this 
phase is by dry mergers only, seen as sudden upwards turns 
after Mstar (z) has flattened (Fig. in contrast with the 
smoother growth of M atar due to star formation. 

The most massive galaxies have final stellar masses in 
the range 10 n ' 7 M Q < A/star £ 10 11,9 M S . In these galaxies, 
M s tar 3> Mg t r a * because their haloes started forming very 
early and crossed M cr i t when the entry mass was larger 
than M s c t r l r ~ 10 n M Q . The upturn of Merit (z) at z > 3 
causes most of these giant galaxies to become red-and-dead 
almost simultaneously at z ~ 3.5 — 4. Their post- quenching 
luminosity evolution is reasonably well described by a simple 
relation of the form 

M B 21.5 -z ±0.5 (4) 

(diagonal solid line and ±0.5 mag dotted lines in Fig. |4]). 
The haloes of these galaxies correspond to the sites of galaxy 
clusters at low redshifts. 

Frequent mergers inside massive haloes provide an ad- 
ditional mechanism for the growth of giant galaxies that 
can raise their mass up to M s tar ~ 10 12 Mq. Indeed, five 
out of our six massive galaxies grow by a factor of ~ 2-3 
through dry mergers after entering the red sequence. This 
lat e merging agree s with semi analytic modelling ca l culati ons 
bv lDe Lucia et all |200rj) and |Pe Lucia fc Blaizotl (|2007f ). A 
note of caution regarding the comparison of these predic- 
tions to observations is that the light associated with the 
growing stellar mass at the centres of galaxy clusters may 
be spread out across a non-negligible fraction of the clus- 
ter volume. Combined with observational limitations, this 
may lead to an underestimate of the stellar mass in bright- 
est cluster galaxies at the centers of cluster s (|Gonzalez et all 
l2005t lLauer et al]|2007l ; iFaber et al]|2007f ). 

The upper-intermediate mass bin in Fig.|4](10 11 ' 2 MQ < 
M s tar < 10 3 Mq) is populated by central galaxies in which 
Mhaio has crossed M cr i t at z 3, where M cr it = M s hock 
(Eq. 1). All of these galaxies enter the red sequence with 
a same stellar mass of Mg* ~ 10 n M Q . In both the high 
and the upper-intermediate mass bin there is a strong cor- 
relation between the final stellar mass and the redshift at 
which the galaxies become red. Galaxies that enter the red 
sequence earlier end up being more massive at z ~ 0, even 
when the entry mass is similar, because they have more time 
to grow through dry mergers along the red sequence. 

In contrast, the two lower mass bins (1O 1O ' 7 M0 < 
Mstar < lO lo - 8 Af and 10 10 - 2 Af o < M Btar < 10 103 M Q ) are 
dominated by satellite galaxies, identified in Fig. [4] by an 
abrupt big rise in their halo mass growth as they fall into a 
massive host halo. As satellites, they have a lower merging 
rate, with the growth by dry mergers along the red sequence 
playing a role only in a couple of the low-intermediate mass 
galaxies, and in none of the low-mass galaxies. For most of 
the low-mass galaxies, the transition to the red sequence 
corresponds to a merger of their dark matter halo into a 
halo more massive than M cr it, causing them to shut down 
as satellites. As the timings of these events are not strongly 
correlated with the halo mass before merging, the low-mass 
galaxies display a wide spread in quenching times, which are 
not directly related to final mass. 
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Figure 5. Time evolution in the colour- mass diagram for galaxies that are currently in four different mass bins along the red sequence, 
which are the same as in Fig.|3]and Fig. [4] The history tracks refer to the most massive progenitors. The colour assigned to each galaxy 
is the same as in Fig. [4] The green dashed line marks the green valley that separates the red sequence from the blue cloud at z ~ 1 in 
our model (see Fig- [6j) . The symbols refer to two different redshifts, as indicated, and show the pace at which galaxies move along the 
colour-mass tracks. Ignoring bursts of star formation, galaxies get steadily more massive and redder until they turn much redder in a 
short time once they are quenched. Small satellites are quenched and then keep a fixed mass. Massive galaxies continue to grow by dry 
mergers after they are quenched. The massive galaxies with Mstar Si 10 117 Mq have already been quenched by z = 3. Dry merging is 
not apparent at Mstar ~ 10 1:l Mq. It becomes an important growth mechanism above this mass. We have plotted this figure assuming 
that all galactic discs are observed face-on. 
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One lower-intermediate mass galaxy (light blue) resides 
in a halo less massive than M cr i t (and thus appears as one 
of the 7% blue symbols in Fig. [3} . This galaxy has not been 
quenched by any of the three explicit mechanisms, but rather 
ran out of gas by continuous star formation. 

One of the low-mass galaxies (fuchsia) has been the cen- 
tral galaxy of a halo that reached nearly 10 13 Mq before it 
crossed M cr it at z ^ 3. However, this galaxy remained small 
because it was quenched before it had time to convert into 
stars the gas that had accreted and was accreting onto it. 

To conclude, Fig. [4] demonstrates that shutdown at a 
constant critical halo mass introduces a rather constant en- 
try mass for central galaxies into the red sequence. Central 
galaxies that enter the red sequence earlier end up more mas- 
sive because they have more time for subsequent growth by 
dry merging. Not e that the downsizing of M cr j t at high red- 
shift predicted bv lDekel fc Birnboiml ( |2006l ) in itself induces 
a downsizing in the entry mass, which enhances the effect 
of archaeological downsizing in today's giant ellipticals. Dry 
merging is important in the growth of massive red galax- 
ies with M s tar 10 Mq. Galaxies that become satellites 
enter the red sequence with masses lower than the charac- 
teristic entry mass, and they dominate today's population 
at M s tar < lO n M . 



4 DOWNSIZING IN THE ASSEMBLY OF THE 
RED SEQUENCE 

In this section, we show that the shutdown above a critical 
halo mass produces downsizing in the buildup of the red se- 
quence, in the sense that the high-mass end is populated at 
earlier times. For this purpose, we consider in Fig. [6] a large 
randomly selected sample of 8,000 galaxies, except that all 
the giant galaxies with M ata r > 1O 11,5 M0 are included for 
better statistics. We plot a series of colour-mass diagrams, 
which show these galaxies at z = 0, and their main progeni- 
tors at higher redshifts. The red symbols highlight the galax- 
ies that end up on the red sequence in four distinct mass in- 
tervals, LogM star = [10,10.5], [10.5,11], [11,11.5], [11.5,12]. 

The overall population exhibits the same behaviour seen 
in Fig. U The distinction between a blue cloud and a red se- 
quence is alre ady present at z ^ 2 (i n agreement with FIRES 
observations; [Giallongo et al.ll20"05T ). with about 60% of the 
most massive galaxies already on the red sequence at that 
redshift. Upper intermediate mass galaxies grow along the 
blue cloud until they reach the entry mass, near the upper 
mass limit for the blue cloud, where they turn red. The most 
massive galaxies climb up to the top of the blue sequence, 
turn red, then continue their growth by dry merging along 
the red sequence. 

The mass growth due to dry mergers is a factor of ~ 2.5 
in the highest mass bin and ~ 1.7 in the upper- intermediate 
mass bin. It becomes negligible below 1O 11 M0. This tran- 
sition from wet to dry mergers may be related to the ob- 
served division of elliptical galaxies into two major types, 
giant boxy non-rotating ellipticals versus the smaller discy 
rotating ellipticals. Indeed, the origin of this distinction has 
been proposed to be the difference between dry merging 
along the r ed sequence and r e cent wet mergers along the 
blue cloud (IBender et al.lll988l. fl989l; iNieto fc Benderll 19891; 



iNieto et al.lll99ll ; iBender et al.lll992l ; iKormendv fc Benderl 



ll996tlTremblav fc Merrittll 199rj ; [Gebhardt et al1ll996h . The 

transition from wet to dry mergers may also be related to 
the division of ellipticals according to their inner density pro- 
files, where giant ellipticals show cores and less massive ones 
show power-law cusps, possibly refle cting gas-rich mergers 
(|Faber et al.ll 19971 ; lLauer et alj|2007h . This transition is ob- 
served at M v ~ -20.5, which is ~ lO ia8 M (assuming 
M/L v ~ 6 x lQ- - 092(Mv+22) M(7)/ La based on the M/L es- 
timates of iGebhardt et ail | |2003l )). This is indeed close to 
our estimated entry mass of M^^l ~ 10 1 M , below which 
the growth is mainly by wet mergers. 

The same transition mass is typically the minimum en- 
try mass for central galaxies. It is also comparable to the 
bright edge of the blue cloud. The fact that this bright edge 
is constant with time reflects the fact that the entry mass 
for central galaxies is roughly constant after z ~ 2. 

The low-mass bin and the lower intermediate mass bins, 
which are dominated by satellites, are mainly populated af- 
ter z ~ 1. Low mass galaxies make an abrupt transition 
from blue to red at masses lower than the minimum entry 
mass for central galaxies. This is consistent with the notion 
that they are quenched when they fall into larger haloes and 
become members of groups or clusters. 

We notice that at z ~ 2 the red sequence has a rather 
constant colour independent of mass, and that a tilt develops 
gradually in time. The origin of this effect, and its relation 
to downsizing, will be addressed elsewhere. 



5 ARCHAEOLOGICAL DOWNSIZING IN 
STELLAR AGES 

Perhaps the most useful way to measure downsizing is by 
comparing the distribution of s tellar ages in galaxies of dif- 
ferent masses. iThomas et ail (|2005l ) deduced stellar ages 
from the absorption line indices H/3, Mg b and <Fe> for a 
sample of 124 early-type red sequence galaxies. They fitted 
the distribution of stellar ages in galaxies of a given stellar 
mass using a Gaussian, with the mean time and standard 
deviation as free parameters. Their main finding, as read 
from their Fig. 10, is that the mean stellar age is increas- 
ing with mass, while the scatter is decreasing with mass. 
They also find that the mean stellar age is increasing with 
environment density at a given stellar mass. It is important 
to realize that the measured age spreads are inferred from 
[a/Fe] and not from detailed modelling of galaxy spectra. 
Moreover, the spread in the times at which different galaxies 
attain their peak SFR may be significantly larger than the 
age spread in each galaxy, which is what [a/Fe] measures. 
Taking these differences into account is important when we 
compare the age spread estimated by I Thomas et all (|2005l ) 
with those obtained from the mean star formation histories 

in our model. 

T o compare our model with the results bv lThomas et all 
(|2005l ). we should consider all the stars that end up in a 
given red galaxy, and not just those that formed in its most 
massive progenitor. Fig.[7]plots the predictions of our model 
for the distribution of stellar ages in the spheroids of bulge- 
dominated galaxies. We compare the results for the standard 
model (no shutdown, except in satellites) and the new model 
(shutdown of star formation in satellites, in massive haloes 
and in bulge-dominated galaxies). Once again, we split the 
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Figure 6. Evolution of the colour-mass diagram from z ~ 2 to z ~ 0. The grey dots represent 8,000 randomly selected galaxies at z = 
(and all the giant galaxies), and their main progenitors at higher redshifts. The galaxies with M s t ar > 10 10 Mq on the red sequence at 
2 = are divided into four bins according to their final stellar mass, as marked by the red symbols in the four columns of the upper row. 
The red symbols in the lower panels highlight the corresponding main progenitors at the different redshifts. 
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Figure 7. The distribution of stellar ages as a function of mass for bulge-dominated galaxies, obtained by averaging the specific SFR 
histories in the given mass bins. The specific SFR is inferred from the age distribution of all the stars in today's galaxies (not just those 
in their main progenitors). Left: the standard model. Right: the new model, in which star formation is shut down above a critical mass 
and in bulge-dominated central galaxies. 



red galaxies into four stellar mass bins. We see that, in 
the standard model, galaxies of all masses have similar star 
formation histories. The weak downsizing signature proba- 
bly reflects the natural downsizing of dark haloes in which 
star formation occu rs only above a minimum halo mass 
l|Neistein et al.ll2006t ). In contrast, the new model shows a 
strong downsizing effect, where more massive galaxies form 
stars earlier and over a shorter time span, as observed. 

In both models, the SFR histories for all masses are 
predicted to have long tails extending to the present time. 
This has also b e en see n in the semi-analytic simulations by 
lDe~ Lucia et al.l (|2006t ). Note, however, that the curves in 
Fig. [7] are averages over many galaxies. Individual galaxies 
have a narrower age spread. The general skewness of the 
SFR toward early times reflects the overall trend of the star 
formation history in the Universe, where roughly half the 
stars are made before a look-back time of t = 8 Gyr [z = 1). 
In the new model, the SFR in massive galaxies is truncated 
by crossing M cr i t early. Less massive central galaxies pass 
through Merit later, so their star formation stops later. 

The two lowest mass bins have the widest age spreads 
because satellites fall into bigger haloes and get quenched 
over a wide range of times. Small differences between the 
two lowest mass bins reflect both the residual effect of central 
galaxies in the lower intermediate mass bin and the tendency 
of more massive satellites to fall into massive haloes earlier. 

Fig. [8] compares the results of the new model with 
th e mean stellar ages and internal age spreads inferred 
bv lThomas et all (|2005l ). including the environment depen- 
dence. The model largely reproduces the observed effects. 
In the two environments, the model mean stellar ages in- 
crease with galaxy mass in general agreement with the ob- 
served trends, except that the model has a slight tendency 
to overestimate the stellar ages in the field (see below), but 
the big difference between environments may be an artifact 
of the small sample used (D. Thomas, private communica- 
tion). In our model, the stellar populations of the most mas- 
sive ellipticals are as old or even slightly older than those in 
rTh et al.l (|2005h . even in the densest environments. The 
model has no trouble making old-enough stars in massive 



15 






10.0 



GallCS Thomas et al. 05 . 

M,.»,„<M„ rf , __. Field 

_ Cluster 



"halo^ m crlt 

M hn i_>M_ rI . 



10.5 11.0 11.5 

Lo g( M bulge/ M ©) 



12.0 



Figure 8. The median stellar age of bulges in bulge-dominated 
(E/S0) galaxies as a function of mass. The symbols and error bars 
correspond to the look-back time by which (50±34)% of the bulge 
stellar mass has formed. To mimic environment density, the green 
and black symbols refer to galaxies in haloes below and above the 
critical mass ~ 2 X 10 12 Mq, respectively. The red lines show in 
comparison the ±<r age spread as inferred bv lThomas et alj 1120051 ) 
for E/S0 galaxies in field and cluster environments. 



galaxies. At any given galaxy mass, the age of the stellar 
population is higher in a cluster environment than in the 
field. 

The internal spread in stellar ages is narrower for more 
massive galaxies, both in the model and the data. While in 
the model the most massive galaxies formed all their stars in 
a short time spa n 2 Gyr, thi s time span is not as brief as 
that inferred by iThomas et "all (|2005l) from their measured 
values of [a/Fe], which is ~ 0.3 Gyr. At a time of t ~ 2 Gyr 
after the Big Bang, the characteristic dynamical time for 
flowing at the virial velocity from the virial radius to the 
halo centre is ~ 0.4 Gyr. The observed duration thus implies 
a massive burst of star formation on the shortest possible 
dynamical time while the model predicts a longer duration, 
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which is comparable to the Hubble time at that epoch. A 
burst as short as inferred by iThomas et"aH (|2005t ) requires 
a process that accelerates the SFR and then rapidly shuts 
it down, as it could possibly be obtained by major merg- 
ers, q uasar feedback, or shocked accretion |Birnboim et al.l 
l2007h . Such a process is not incorporated in our current 
model. 

The observed mean stellar ages are expected to be 
slightly lower than the model predictions for several rea- 
sons. First, the observed light- weighted average ages are 
younger than the model mass- weighted averages because the 
younger stars are brighter and dominate the light. Further- 
more, younger populations have stronger Balmer equivalent 
widths, which dominate the total equivalent widths, and 
thus add to the brightness effect just mentioned, making 
light-weighted ages still younger. Finally, young tails in the 
age distribution, if present in real galaxies, would pull the 
peak of the Gaussian fit toward younger ages. Considering 
these limitations, the agreement between model and data is 
good. 



6 DOWNSIZING IN STAR FORMATION VS. 
UPSIZING IN ASSEMBLY 

In Fig. [9] we extend the analysis of the new model to other 
times that characterize the SFR and its shutdown. These 
are the times (a) of peak SFR, (b) of shutdown, (c) at which 
Afhaio crosses M cr i t , (d) at which the galaxy becomes bulge 
dominated, and (e) at which the main progenitor attains 
half the final galaxy stellar mass. Each characteristic time 
can be interpreted as showing either downsizing or upsiz- 
ing depending on whether it is monotonically decreasing or 
increasing as a function of final stellar mass. We shall see 
that all the above time indicators referring to star forma- 
tion show a general downsizing trend, as opposed to the up- 
sizing trend that the model predicts for the analogous time 
indicators that characterize mass assembly. 

It is worth noting the differences between Fig. [9] and 
Fig. [8] First, the times t in Fig. [9] are from the Big Bang, 
while in Fig.[8]they were lookback times, that is to—t, where 
to = 13.7 Gyr is the present age of the Universe in the cos- 
mology used for this study. Second, Fig. [8] is constructed 
from a small number of galaxies, where the error bars cor- 
respond to the internal stellar age spreads in the individual 
galaxies, while Fig.[9]shows the entire red galaxy population 
in the new model. Third, the galaxies in Fig. [8] are all red 
bulges of bulge-dominated galaxies, while Fig. [9] shows all 
the red galaxies, and the times contain the contribution of 
the disc component if present. Fourth, Fig. prefers to all the 
stars in today's bulges, while Fig. [9] refers only to the stars 
in the galaxies' main progenitor back in time (the purple 
symbols in Fig. [9] are an exception). 

The time of maximum SFR, t maX sfr, is shown in Fig.(9]a.. 
In most red galaxies, it occurred within the first ~ 4 Gyr 
from the Big Bang. There is a clear but weak downsizing 
trend in the epoch of maximum SFR. In the most massive 
galaxies, the SFR peaked at t ~ 1 — 2 Gyr, while at low 
masses there is a larger spread in t maxs f r . We thus recover 
downsizing features similar to the 'archaeological' downsiz- 
ing obtained from the stellar age distribution in red galaxies 
at z ~ 0. 
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Figure 9. Times characterizing the SFR and its shutdown versus 
final stellar mass, (a) The time of peak SFR i maxs f r . (b) The 
time of star formation shutdown t sriU tdown ■ ( c ) The time at which 
A4alo crosses M crit , *Mhalo>Mcrit- ( d ) Thc timc at which the 
bulge mass becomes larger than the disc mass tb u i ge >disc- ( e ) The 
assembly time at which the main progenitor attains 50% of the 
final stellar mass, t 5 Q% (red symbols). Also shown is the time at 
which half of the stars have formed (purple symbols). The times, 
measured from the Big Bang, are evaluated by tracking the main 
progenitor of each galaxy back in time. Central galaxies today are 
marked green and satellites arc marked black. Red symbols with 

error hars are median rallies and lower and nnnor rmartiles in ten 
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Figure 10. Role of the two quenching criteria. Each point shows 
the stellar masses of a galaxy when it became bulge dominated 
(j/-axis) and when its halo crossed M cr it (a;-axis). Galaxies that 
crossed M cr n before (after) they became bulge-dominated lie 
above (below) the diagonal red line. Green and black points refer 
to central and satellite galaxies at z = 0, respectively. Galax- 
ies that never become bulge-dominated, mainly satellites, are put 
near the top edge of the plot. Galaxies that never cross the criti- 
cal halo mass, mainly centrals, are put near the right edge of the 
plot. 

For a fixed final stellar mass, the SFR peak occurs ear- 
lier in satellite galaxies than in central galaxies. This is be- 
cause low-mass central galaxies are in low mass haloes in 
which star formation does not start early, it is slowed down 
by supernova feedback, and it is not shut down early. Satel- 
lites, on the other hand, much like high stellar mass galaxies, 
can be part of massive haloes already at early times, so they 
can form stars early and also shut down early (SPSS dat a, 
iBlanton et al.| [2006l; our new model. [Cattaneo et al.|[2006t ). 

The time of shutdown in the most massive progenitor, 
^shutdown, is shown in Fig. [5p. In ~ 7% of the red galaxies 
with M s tar > 1O 1O M0 the formation of stars has never been 
quenched by any of the explicit criteria applied in the new 
model. These galaxies have turned red simply because they 
have exhausted their gas reservoir. In addition, ~ 2% of 
the red sequence at M star > 1O 1O M0 is made of starbursts, 
which are red because they are dusty. We have removed 
these unquenched and dusty galaxies from Fig. [9p and from 
the shown median. The general downsizing trend in Fig. |9p 
above ~ 10 10,7 Mq is similar to Fig. [9^,, especially for the 
central galaxies, but the magnitude of the effect is stronger 
for tshutdown- This is because intermediate mass galaxies not 
only reach their peak SFR somewhat later than more mas- 
sive galaxies, but they also have tails of SFR extending to 
late times, z 1 (Fig. [7Jl. The most massive galaxies un- 
dergo an intense burst of star formation when the Universe 
is only ~ 1.5 Gyr old and then they are almost immediately 
shut down. In red galaxies of intermediate masses, the peak 
SFR is at t ~ 3 Gyr and the formation of stars drags on 
for ~ 6 Gyr after the peak. The shutdown time does not 
show downsizing at low masses, where the red population is 
dominated by satellites. 

We recall that our model contains two quenching mech- 
anisms. One is due to the shock heating in haloes once they 
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Figure 11. Growth by dry merging along the red sequence versus 
stellar mass. It is measured by the ratio between the final stellar 
mass of a red galaxy and the mass it had when it was quenched by 
one of the three shutdown mechanisms of our model. Green and 
black points refer to central and satellite galaxies at z = 0, respec- 
tively. The red symbols with error bars show the median values 
of M s tar /^st^T an< ^ * ne l° wer an d upper quartiles in logarithmic 
bin of 0.2 dex in M s t ar . 

become more massive than M cr it , and the other is a result of 
wet major mergers. Mergers are assumed to simultaneously 
produce dominant stellar spheroids and trigger quasars at 
their centres, which then help removing th e gas available 
for sta r formation in merger remnants (e.g., ISpringel et alj 
l2005al ). The bulge-dominance criterion applied in the the 
new model is meant to mimic quasar feedback of this sort. 
We assess the relative importance of these two mechanisms 
in the new model by splitting Fig. [9p into two, show- 
ing separately the times for crossing the critical halo mass 
tMhaio>Mcrit (Fig. \9fc) and for the bulge to become dominant 
£bui g c>disc (Fig. [9}i). Galaxies that have never crossed the 
critical halo mass do not appear in Fig. [5J:, and galaxies 
that have never become bulge-dominated do not appear in 
Fig. 031. 

Massive galaxies with M s tar > 10 11 Mq and ^shutdown < 
3 Gyr in Fig. |9p populate the same region of the diagram 
in Fig. [9li and Fig. [9JI. These are central galaxies of massive 
haloes, which simultaneously obey Mhalo > Merit and have a 
dominant spheroid due to intense merging. Galaxies that are 
less massive tend to show up either in Fig. [9ji or in Fig. [9}1, 
depending on whether they are satellites or small central red 
galaxies. 

The same three populations can be identified in Fig. 1101 
which displays for every galaxy the stellar masses at the two 
characteristic times: when it entered a halo more massive 
than the critical mass, and when it became bulge-dominated. 
The cloud of black dots at the very top consists of satel- 
lite galaxies in massive haloes, which have never become 
bulge-dominated because the merger rate is low for satel- 
lites. This population corresponds to the black symbols in 
Fig- Op that survive in Fig. [9J; but not in Fig.[9jl. The cloud 
of green dots at the very right consists of bulge-dominated 
central galaxies in haloes that are just below the critical 
mass. These galaxies correspond to the intermediate mass 
galaxies shown as green symbols in Fig. [3] Without the 
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bulge shutdown criterion, these galaxies would be on the 
blue cloud, and they w ould extend it to higher luminosities 
(|Cattaneo et alj|2006l ). This population corresponds to the 
green symbols in Fig. [9)3 that survive in Fig. (9ji but not in 
Fig. [5J;. The third population comprises galaxies that both 
inhabit massive haloes and are bulge-dominated. This popu- 
lation contains the most massive objects and is composed of 
both central and satellite galaxies. The distinction between 
central and satellite galaxies is not very useful here because 
many satellites were central until recently. A more useful dis- 
tinction is between the galaxies that were first quenched by 
the bulge criterion and the galaxies that were first quenched 
by the halo criterion (points below and above the diagonal 
respectively). These galaxies obey the two quenching crite- 
ria almost back-to-back because the two events are related 
to the hierarchical growth of dark matter haloes. 

The characteristic times considered so far are indicators 
of star-formation history. Fig. [9^ displays instead an indica- 
tor of the assembly history, the time t 50 % mass at which 50% 
of the final stellar mass has been assembled in a single ob- 
jects (red symbols). The assembly time is roughly constant 
with mass for low and intermediate-mass galaxies, and it ac- 
tually reveals upsizing at the high-mass end. This resembles 
the asse mbly upsizing of the m ain progenitor of dark matter 
haloes in lNeistein et all (|200rJ) . 

In galaxies of low or intermediate mass, t 50 % maaa is com- 
parable to t max sf r of Fig. The complete shutdown of star 
formation occurs when the Universe is ^8 Gyr old, but 
^ 80% of the final stellar mass is already in place in the 
most massive progenitor after the first ~ 5 Gyr. Only little 
mass is added to galaxies of M sta r < 10 1 Mq after this point, 
as the merger rate is not high in this mass range, as demon- 
strated by Fig. 1111 which shows the mass growth factor after 
entering the red sequence as a function of the final stellar 
mass. Moreover, accretion onto the haloes of small galaxies, 
and therefore onto small galaxies themselves, is suppressed 
by tidal effects in the vicinity of more massive haloes (Hahn, 
Dekel et al., in preparation). 

In contrast to small galaxies, galaxies with 
M s t a r 10 1 M© experience substantial growth through dry 
mergers after they cease to form stars (Fig. 1 11 p . Therefore, 
the most massive galaxies continue gathering mass for a 
much longer time. When the Universe is ~ 7 Gyr old, they 
have barely put together half of their final stellar mass. 

The onset of upsizing in assembly is right above the 
minimum entry mass for central galaxies, near the mass 
separating the boxy and discy ellipticals (Section 4). Our 
model hints to a simple explanation for this bimodality of 
structural and kinematic properties of galaxies along the 
red sequence. Discy ellipticals are the products of recent 
wet mergers of gaseous disks along the blue cloud, in haloes 
of Mhaio < Merit. Boxy ellipticals, on the other hand, re- 
sult from dry mergers along the red sequence, typically at 
the centres of haloes with Mhaio > M cr it- The most mas- 
sive ellipticals, which have only finished assembling 80% of 
their final stellar mass in the last few billion years, are to be 
identified with the brighte st cluster galaxies a t the c entres 
of Abell clusters. Also see Ide Lucia fc Blaizotl i|2007t ) for a 
semi-analytic study. 

Fig. [5£ also displays for comparison the time at which 
50% of the final stellar mass was formed (purple symbols), 
showing a star-formation downsizing similar to the trends 
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Figure 12. The evolution of the mass functions of red and blue 
galaxies from z = 1 (dashed lines) to z = (solid lines) in our 
model with star formation shutdown. 



seen in the other panels. This figure demonstrates how up- 
sizing in assembly coexists with downsizing in star formation 
in the same hierarchical clustering scenario. 

A second final remark on Fig. [9] concerns the difference 
between central and satellite galaxies. Compared to satel- 
lites of the same mass, in central galaxies the SFR peaks 
later, it shuts down later, and the stellar mass assembles 
later (compare the green and black dots in Fig. [9^, b and 
e). 



7 DOWNSIZING IN THE GALAXY MASS 
FUNCTIONS? 

The time evolution of the mass functions of red and blue 
galaxies has been used as another indicator for downsizing 
in the formation of the red sequence. The features consid- 
ered, for example, are (a) the downward shift with time of 
the mass where the red and blue functions cross each other, 
and (b) the indications for a late increase in the mass func- 
tion of red galaxies at low masses compared to an early in- 
crease at high masse s (|Bundv et al 1 l2005l ; lBorch et al.ll2006r 
Pannella et ahlbOOrj ; also see the discussion in iFaber et al.l 
2007h . In this section we highlight uncertainties in such argu- 
ments, both due to inherent difficulties in interpreting model 
mass functions and to the large uncertainties that still exist 
in the data themselves. Our principal conclusion is that the 
mass functions of red and blue galaxies are poor indicators 
of downsizing. 

The observational uncertainties are discussed separately 
in the Appendix. Here we concentrate on the predictions 
of the model and demonstrate how hard it is to use mass 
functions to infer downsizing or upsizing, even if the data 
were perfectly accurate. 

The evolution of the mass functions of red and blue 
galaxies in the new model from z — 1 to z — is shown in 
Fig. [12] The red mass function has grown at all masses, while 
the blue function has decreased or remained constant. Below 



by ~ 0.2 dex at all masses, accompanied by a comparable 
uniform decrease in the blue function over the same mass 
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range. Above M s c t "£, the red function grows by a larger factor 
while the blue function remains roughly constant. 

Our earlier analysis can help us identify the processes 
responsible for the buildup of the red sequence at the low- 
and high-mass ends. Red galaxies with Mstar > Mg t " r consist 
of stars that have formed at z 2> 1, so their progenitors were 
already on the red sequence at z ~ 1. These galaxies have 
continued assembling through dry mergers of red galaxies, 
such that the growth of the red population at the high-mass 
end comes largely at the expense of the red population at the 
low-mass end, from below M s c t "'. The true influx of galaxies 
into the low-mass red population equals the increase of the 
red mass function seen in Fig. [T2] plus the transfer of red 
galaxies to the high-mass end by dry mergers. This influx 
must come from blue galaxies migrating to the red sequence 
because the increase in total stellar mass in red galaxies 
cannot come from new star formation in the red sequence 
jFaber et alj|2007h . 

Part of the stellar mass transferred from the blue to the 
red sequence is readily accounted for by the decrease in the 
number of blue galaxies at Mstar < M^axj where most of the 
blue sequence mass resides. However, the total stellar mass 
in galaxies with M s tar > 1O 1O M0 has increased from z ~ 1 
to z ~ 0. This implies that, despite the drop in the number 
of blue galaxies, new stellar mass must have been formed 
along the blue sequence at £ z £ 1. In our model, this is 
mainly due to the conversion of disc gas into stars. Only to 
a lesser extent is it due to new galaxies moving into the blue 
sequence due to growth through the accretion of gas and 
wet mergers. 

We have thus identified three fluxes that drive the evo- 
lution of the populations of red and blue galaxies in the 
colour-mass diagram from z ~ 1 to z ~ 0. First, the flux 
of stellar mass into the blue sequence due to conversion of 
gas into stars in blue galaxies £ M%£- Second, the tran- 
sition of blue galaxies onto the red sequence at £j M^f„. 
Third, the transfer of red galaxies by dry mergers from the 
low-mass part of the red sequence to the high-mass part, 
M s tar > Mt$£. Only the first mode generates new stellar 
mass in the overall galaxy population, and only the first two 
contribute to the growth of stellar mass in the red sequence. 

Now, with the benefit of hindsight, we are ready to in- 
terpret the changes seen in Fig. 1121 One obvious trend in 
both the model and the figure is the assembly upsizing that 
creates the most massive elliptical galaxies on the red se- 
quence at late times. This is visible in figure as the rightward 
shift in the red function at high masses. These galaxies have 
assembled late and therefore appear at high masses late. 
Less clear, however, is the behaviour of the red mass func- 
tion below M fi c t j*. We know from the detailed histories of 
galaxies in Fig. [9] that the star-formation downsizing trend 
at these masses is strong, yet this effect is not obvious in 
the mass functions. All masses below M^J appear to grow 
about the same. To be precise, there does not appear to 
be much greater growth at the lowest masses on the red 
sequence compared to the growth at M^ T . This happens 
largely because our model contain many satellite disc galax- 
ies that form too early and are quenched too early (Fig. IA1I 
in the Appendix). The lesson is that deducing either down- 
or upsizing trends from the mass function is tricky. 

Fig. [12] also sheds light on a common misinterpreta- 
tion of the crossover mass, that is, the point where the red 



function intersects the blue function. Downward movement 
of this point in time has been s uggested as an indic ator of 
downsizing (|Bundv et all 120061 : iHopkins et aT]|2007h . Even 
though downsizing is actually taking place at these masses 
in the model, we have seen that this is not evident from each 
of the mass functions alone. Instead, Fig. [T2] shows quite 
clearly that the crossover point is moving downward sim- 
ply because the number of red galaxies is growing relative 
to the number of blue galaxies, independent of how these 
changes depend on mass. This effect is enhanced when the 
mass functions are relatively flat, which is true both in the 
model and in real data (see the Appendix). The crossover 
point would move downward even when blue galaxies be- 
come red in a uniform rate at all masses, namely without 
any real downsizing trend in the SFR along the blue cloud 
or in the buildup of the red sequence. 

In the Appendix, we make three points concerning the 
observed mass functions. First, determining the number of 
very massive galaxies at the steep bright end is tricky be- 
cause even small random mass errors create large errors in 
the counts. The observed mass function at the bright end is 
not the true function but rather the true function smoothed 
by errors. Owing to the steepness of the bright end, this 
tends to falsely inflate the numbe r of massive objects (the so- 
call ed Eddington effec t; see also iKitzbichler fc White! |2006| 
and lLauer et al.l 1200 7l ) . The Appendix shows that random 
rms errors of only 0.2 dex in mass inflate the observed counts 
by a full order of magnitude at M sta r ~ 10 M§. Model 
counts must therefore be smoothed to mimic this effect be- 
fore comparing to data, but it is unclear whether the errors 
get smaller or larger with redshift, and thus whether the 
evolutionary effect is positive or negative. 

Second, the total magnitudes of massive ellipticals may 
be systematically underestimated because the light profiles 
of very massive galaxies are nearly isothermal and it is hard 
to estimate the total light, let alone distinguish galaxy light 
from intracluster light. This problem likely worsens at low 
redshifts, and so the sense is probably to underestimate 
galaxy masses, and thus the number of massive galaxies to- 
wards low redshifts. 

Third, before comparing mass functions from different 
authors, it is necessary to convert them to the same assumed 
stellar initial mass function (IMF) . Certain authors who for- 
merly appeared to agree now disagree after this correction 
is made. Consistent trends in even just the raw counts are 
not yet evident. 

Given these difficulties, both observational and theo- 
retical, we do not attempt in this paper to deduce either 
down- or upsizing trends from luminosity functions or mass 
functions. All that can be said is that the number of red 
galaxies has increased by at least a factor of 2 near L* and 
M, since z = 1, while the number of blu e galaxies has not 
increased since z = 1 |Faber et al.ll200?f ). These trends do 
not alone argue for down- or upsizing, but at least they are 
in reasonable agreement with the results of our model. 



8 DISCUSSION AND CONCLUSION 

We showed that the downsizing of elliptical galaxies, where 
the more massive galaxies formed their stars earlier and over 
a shorter period, is not in conflict with the standard hi- 
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erarchical clustering scenario. We demonstrated that this 
is in fact a natural outcome of the shutdown of star for- 
mation in haloes above a critical mass Merit ~ 10 12 M Q . 
Such a shutdown is essential for explaini ng the distribu- 
tion of galaxies in lum i nosity and colour (ICattaneo et al.l 
l200d : iBower et ai]|2006l ; ICroton et all 120061 ). In particular, 
this shutdown prevents the overgrowth of the central galax- 
ies in massive haloes beyond the observed upper limit for 
bright blue galaxies, and allows the appearance of mas- 
sive red and dead galaxies, thus reproducing the observed 
galaxy bimodality. The critical mass for shutdown originates 
from stable virial shock heating, which occurs in haloes 
above this threshold and shuts down the cold ga s supply 
for s t ar formation ther e (Birnboim & Dckcl 2003|; iBinnevI 



|2004 iKeres et alj|2005l ; iDekel fc Birnboiml 2006). This see 



nario only works if gas in massive haloes is maintained 
hot by a continuous energy source that balances the radia- 
tive cooling over long cosmological periods. The mechani- 



ing candidate ((Fabian et al.l 2003: 


Ruszkowski et alj 


2004; 


Forman et alj|2005: Best et alj|2005 


;IVoit & Donahue! 


2005; 


Best et al. 20061: iDunn & Fabian 2006). which can result in 



self-regulated accretion and AGN activity capable of main- 
taining the required long-term shutdown (lOmma fc Binnevl 
|2004 iRaffertv et"aU 120061 ; ICattaneo fc Tevssierl 120071 ). An 
alternative/additional source of long-term quenching is the 
gravitational energy of accretion into the centres of the po- 
tential wells of dark matter haloes llBirnboim et al. I l2007t 
IDekel fc Birnboimll2007l ; iKhochfar fc Ostrikerll2007l ). 

The shutdown at a critical halo mass introduces an en- 
try mass Af,rtar to the red sequence. Central galaxies cease 
to make stars and turn red and dead at a mass equal to 
the characteristic mass of the central galaxy in a halo of 
mass M cr it ■ The central galaxies of the most massive haloes 
reach M££ ~ 10 n M© earlier (Fig.gJ. Therefore, they have 
more time to grow above M s c t "r by dry merging along the red 
sequence. For this reason, the central galaxies of the most 
massive haloes are the most massive galaxies in the Universe 
and those that contain the oldest stellar populations. In our 
model, the early growth of the central galaxies of the most 
massive haloes is also helped by the fact that M cr it is as- 
sumed to be high er than 2 x 10 12 Mp before z Si 3, based on 
the predictions of IDekel fc Birnboiml (|2006l ). 

We demonstrated that shutdown-driven downsizing is 
in agreement with the observed archaeological downsizing 
of elliptical galaxies, as measured by Thomas et al. (2005; 
Fig. |8]) . The introduction of shutdown above a threshold 
mass also explains why the emergence of the red sequence 
has started already at z ~ 2 and why the shutdown be- 
gins at high masses (Fig. |6]). Over 3/4 of the model ellip- 
ticals with Mstar > 3 x 10 11 M© were already on the red 
sequence by z ~ 2, in agreement with FIRES observations 
(|Giallongo et alj|2005l) . In contrast, the buildup of the red 
sequence at M sta r 10 11 M Q continues at z ^ 1. 

The downsizing in star formation and the upsizing in 
mass assembly affect the evolution of the mass functions of 
red and blue galaxies. The three phenomena that control 
the evolution of the mass functions are: the formation of 
new stellar mass in blue sequence galaxies, the migration of 
blue galaxies onto the red sequence, and the migration of 
red galaxies from the low-mass end to the high-mass end of 
the red sequence due to dry merging. However, the combi- 



nation of downsizing in star formation and upsizing in mass 
assembly is such that the evolution of the mass functions is a 
poor indicator of downsizing even if the data were perfectly 
accurate (the large uncertainties that still exist in the data 
themselves have been discussed separately in Appendix A). 

Giant ellipticals grow by ~ 0.4 — 0.5 dex after entering 
the red sequence and by ~ 0.2 — 0.3 dex between z ~ 1 and 
z ~ 0. Within the red population, they are those that have 
finished to as semble their mass mo st r ecently (Fig. red 
symb ols; also lDe Lucia et alj|2006l and lDe Lucia fc Blaizot] 
l2007l ). despite the fact that their final mass was turned into 
stars earlier than in other galaxies (Fig.[U£, purple symbols), 
due to the early shutdown of star formation in these objects. 

In this paper we have focussed on the contribution of 
the shutdown above a critical halo mass to the downsizing. 
In principle, there could be other reasons for the observed 
downsizing. Nevertheless, the comparison of the stellar ages 
in the two models, which differ only by imposing an explicit 
shutdown, provides clean evidence that a shutdown of this 
sort can be the main driver of the observed downsizing. 
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APPENDIX A: GALAXY MASS FUNCTIONS 
ATX~0 AND Z ~ 1 

The evolution of the mass functions of red and blue galax- 
ies, whether increasing more strongly at low mass or at high 
mass, has been taken as a barometer of whether the red 
sequence forms primarily through downsizing or upsizing 
jBundv et al.ll2005l ; iBorch et al.ll2006l: IPannella et al.|[2006T 
also see the discussion in lFaber et al.ll2007l ). In Section 7, we 
have shown how the evolution of the mass function of red 
galaxies is driven by a combination of star formation down- 
sizing (due to quenching) and stellar mass assembly upsizing 
(due to merging). We have argued that the inherent difficul- 
ties in interpreting model mass functions, together with the 
large uncertainties that are still present in the models them- 
selves, should suggest caution in using mass function data to 
infer up- or downsizing. In this Appendix, we consider the 
observational data on mass functions at z ~ and z ~ 1. 

The SDSS and 2MASS are the main data s ets used to 
study the mass function of galaxies at z ~ (e.g. iBell et all 
120031 ). Meanwhile, surveys that are increasingly deep (MU- 
NICS, FDF, GOODS) are beginning to give us estimates 
for the mass function of g alaxies at z ~ 1 (e.g. iBundv et alj 
l2006l ; |Panneila" et al. 1 120061 ). Observational determinations of 
galaxy mass functions are affected by three sources of error, 
which we have enumerated in Section 7, and which we now 
discuss in greater detail. 

The first is the steep slope of the mass function at the 
bright end. Therefore, even small random mass errors cre- 
ate large errors in the counts. The observed mass function 
at the bright end is not the true function but rather the 
true function smoothed by errors. Owing to the steepness 
of the bright end, this tends to falsely inflate the number 
of massive objects (the so -call ed Eddington effec t; see also 
iKitzbichler fc White! [20061 and lLauer et all 120071 ). Random 
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rms errors of only 0.2 dex in mass inflate observed counts by 
a full order of magnitude at M sta r ~ lO n ' 5 M0 (compare the 
solid and the dashed lines in Fig. IA1[) . Model counts must 
therefore be smoothed to mimic this effect before compar- 
ing to data, but it is unclear whether errors get smaller or 
larger with redshift, and thus whether the evolutionary ef- 
fect is positive or negative. Fig. lAll shows that we can make 
our model agree with one data set or the other simply by 
making different assumptions concerning the amplitude of 
the random measurement errors. 

A second source that affects the photometry of the most 
massive galaxies is their having shallow light profiles that do 
not converge to a total magnitude. The magnitudes of very 
bright ellipticals are hard to measure on account of their ex- 
tended envelopes, which may be growing with time due to 
galaxy-galaxy interactions. It is possible that many of them 
extend outside the photometric aperture used to measure 
their luminosities. Simply fitting such galaxies with de Vau- 
couleurs R 1 ^ profiles is not appropriate because many pro- 
files (at least locally ) are even more extended than that (e.g., 
iGraham et al.l l200ll). Total magnitude errors may therefore 
be larger for nearby galaxies than for distant ones. For ex- 
ample, the total luminosities of the brightest cluster galaxies 
of near by Abell clusters can v ary by many tenths of a mag- 
nitude (|Gonzalez et al.l 20051) depend ing on the profile fit- 
ting method used. Lauer et al.l (|2007f ) find errors of up to a 
magnitude in SDSS luminosities for such galaxies. These are 
very large errors when one considers that a change of only 
0.2 mag at 4L» translates to a change in number density 
of a factor of two in the Schechter function. Stellar masses 
have even larger errors because they depend on assumed 
evolution models for the fade in M/Lb, which introduce an 
additional uncertainty of a few tenths of a magnitude. While 
random errors are more likely to be important at high red- 
shift, the problem deriving from the extended envelops of 
massive ellipticals is likely to worsen at low redshift, since 
a lot of the growth by dry merging occurs between z ~ 1 
and z ~ (Fig. [4}. Therefore, the combination of these two 
sources of error is more likely to underestimate the evolution 
of the number of massive early type galaxies from z ~ 1 to 
2 ~ 0. 

Thirdly, before comparing mass functions from differ- 
ent authors, it is necessary to co nvert them first to the same 
assumed IMF. iBell et all (|2003l ) used a 'diet' Salpeter IMF, 
which co rresponds to a Salpete r IMF truncated at low stellar 
masses llBell fc de Jond l200ll ). To convert the masses ob- 
tained [PfumsUf^^d! |200(]) for a standard Salpeter IMF 
into mass e s tha t can be compared with the results by 
IBell et al.l (|2003l) . one must subtract ~ 0.15 dex from the 
masse s determined by IPannella et al.l (|2006l ) . iBundv et al.l 
(|2006h used a Chabrier IMF. The masses estimated with 
the Chabrier IMF are lower than those estimated with a 
Salpeter IMF by ~ 0.25 dex. T herefore, we must ad d 0.1 dex 
to the masses determined by IBundv et al.l (120061) to com- 
pare their mass function with that bv lBell et all (|2003l ). Af- 
ter performing these corrections, we see that the app arent 
agreement between t he figures in IBundv et al.l (|200d ) and 
IPannella et al.l (|2006l ) is simply due to the fact that they 
used different IMFs. For reasons that we cannot easily inter- 
pr et, with the same IM F, the galaxy mass function inferred 
bv lBundv et~aH (|200rj ) contains many more massive galax- 
ies at z ~ 1 than the mass function inferred by Pannella 



et al. (2006; Fig.rJT|. The data bv lBundv et alj (|2006l ) are 
consistent with no evolution at M Bta r ^ 10 11 Mp) in the red - 
shift interval £ z £ 1. The data by IPannella et al.1 (|2006l ) 
are consistent with growth in mass by 0.2 — 03 dex at 
AW £ IO^Mq 

In conclusion, it appears that the present state of mass 
functions at 2 = and 2 = 1 is sufficiently confused to 
prevent firm conclusions. Problems are particularly severe 
at the high-mass end, where random errors seriously inflate 
counts at a given mass and questions persist about how to 
measure total light and mass. Fortunately, uncertainties are 
much smaller near L». Though the data are not perfect, the 
general trend is for red galaxies to increase in number with 
time near the knee of the mass function, by amounts that 
are in broad agree ment with both luminosity function data 
(|Faber et al.ll2007l ) and with the predictions of our model. 
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Figure Al. Mass functions at z = and z = 1 for early type galaxies (red) and late type galaxies (blue) classified spectrally (left) and 
morphologically (right) The coloured and black points with error bars are the data. The lines show the predictions of our new model 
with star formation shutdown without (solid lines) and with (dashed lines) the inclusion of a 0.2 dex Gaussian random error. Note how 
the predictions of the model are completely altered at the massive end when such smoothing is included. 
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